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g '■ ABSTRACT 

The Chern-Simons Ginzburg-Landau theory for the fractional Quantum Hall effect is stud- 
io ! ied in the presence of a confining potential. We review the bulk properties of the model 

and discuss how the plateau formation emerges without any impurity potential. The effect 
^ I is related to changes, by accumulation of charge, at the edge when the chemical potential 

is changed. Fluctuations about the ground state are examined and an expression is found 
1-^ ■ for the velocity of the massless edge mode in terms of the confining potential. The effect 

C ! of including spin is examined for the case when the system is fully polarized in the bulk. 

In general a spin texture may appear at the edge, and we examine this effect in the case 

of a small spin down component. The low frequency edge modes are examined and a third 

order equation is found for velocities which indicates the presence of three different modes. 

The discussions are illustrated by numerical studies of the ground states, both for the one- 

and two-component cases. 



1 Introduction 

Since the experimental discovery of the quantum Hall effect (QHE) |l], 0, which occurs 
in a two-dimensional electron gas subject to a strong perpendicular magnetic field, several 
approaches to a theoretical description of this system have been developed. As a useful 
supplement to existing microscopic descriptions P], Zhang, Hansson and Kivelson 0] in 
1988 proposed a Chern-Simons Ginzburg-Landau (CSGL) model for the fractional quantum 
Hall effect (FQHE). This model is based on the concept of "statistical transmutation", i.e. 
the fact that in two dimensions, fermions can be described as (charged) bosons carrying 
an odd integer number of (statistical) flux quanta. At the Lagrangian level this is done by 
adding Chern-Simons fields coupled to the bosons. In this way, electrons in an external 
magnetic field are described as bosons in a combined external and statistical magnetic 
field. At special values of the filling fraction the statistical field cancels the external field 
(in the mean field sense) and the system is described as a gas of bosons feeling no net 
magnetic field. These bosons "condense" into a homogeneous ground state. In fact, this 
rather simple effective theory reproduces several of the key features of the FQHE |^, such 
as the quantization of the Hall conductance and the existence of (anti)vortex excitations. 
It is even possible to re-construct Laughlin's electron wave functions from this model by 
considering quantum fluctuations about the (mean field) boson ground state p. 

An interesting aspect of the quantum Hall system is the existence of gapless edge 
modes which has important consequences for the transport properties of the system 
0. Several authors have studied the CSGL model in the presence of an (infinitely) steep 
external confining potential and have shown the presence of gapless edge modes P, P, p!0 



in this model. The description of these modes in the CSGL model has also been related 
Tl| to the chiral Luttinger liquid description of the edge excitations ||T2|, Q . 



In this paper we study a number of aspects of the CSGL model, discussing both its 
ground state and excitations. We focus, in particular, the attention on edge effects and 
follow up the discussion of previous papers on this subject by considering the density profile 
and edge excitations for a smooth confining potential. We also examine the effects on the 
edge of including the electron spin in the description. 

A short presentation of the model in its simplest form is given in section ^ It describes 
the quantum Hall system without spin at Laughlin fillings z/ = l/(2m -|- 1). In section ^, 
we review the ground state properties of the bulk and discuss how the existence of finite- 
energy vortex excitations leads to incompressibility of the ground state and the existence of 
quantum Hall plateaux, two key features of the QH system. The edge of the system in the 
presence of a smooth confining potential is then studied (section ^). In the ground state 
(with uniform bulk density) one finds a one-parameter family of solutions parametrized 
by the total charge at the edge 0. When this charge exceeds an upper critical value, the 
ground state with uniform bulk density becomes unstable due to formation of anti- vortices. 
Similarly, there is a lower critial value of the total charge beyond which the bulk tends 
to rearrange itself by formation of vortices. A numerical study of the ground state within 
these limits is presented. We further perform a systematic analysis of edge excitations and 



show the presence of a inassless mode with dispersion of the expected form uo/q =< E > / B 
0, and with an exphcit expression for the averaged electric field < E >. 

In section ^ we generalize the model to include spin by introducing a two-component 
description with a Zeeman interaction. The ground state is completely polarized in the 
bulk, however, a spin texture at the edge may be energetically favourable if the Zeeman 
energy is not too large. We study the ground state analytically and numerically; there is a 
two-parameter family of solutions parametrized by two conserved charges. The numerical 
analysis provides estimates of the critical Zeeman gap and it shows rotation of the spin 
vector along the edge for smaller values of the Zeeman energy. (The same effect has 
been demonstrated in a recent study by Karlhede, Kivelson, Lejnell and Sondhi [|13| who 
used Hartree-Fock techniques and an effective spin model.) We further examine the edge 
excitations of the spin-polarized system and find an equation for the mode frequencies 
which indicates the presence of three massless modes at the edge. 

Finally, in section |^ we summarize and discuss our results. 

2 The model 

The spinless system is described by the Lagrangian [Q 

C = <j)*{tdt - ao - eAo - ^(x) + ;,)0 - «: |(V - z(a + eA)) <j)f - X \<j)\^ + l^e^-^^a^d^a^. (1) 

The Bose field is coupled both to an external field A = (0, Bx) with B < 0, and to 
the statistical field a. /i is a chemical potential, e = — |e| denotes the electron charge, 
6 = (2m + l)iT for the u = l/(2m + 1) state and the term A|0|^ models the Coulomb 
interaction which is present in the underlying microscopic theory. We choose Aq = 0. 
With the Bose field decomposed in terms of density and phase, 

<P=Vp e*'. (2) 

the field equations are 

^^^^Vp =[ao + V-fi + dt9 + K{A- VOf + 2Ap] ^p ,^. 

9,p + V-j = 

and 

V ^A = -2e{p-p) 

(4) 



9x0-0 + dtOx = -'20 jy 
dyao + dtay = 29 j^, 



where ^ = a -|- eA, p = eB/26, and the current density is 

] = -2Kp{A-Ve). (5) 



The Hainiltonian is 



H 



d X 



fi:|D0|2 + r|0|2 + A|0|^-^| 



(6) 



where T) = V — iA and the constraint (^) is iniphcit. 

In some parts of the paper it will be convenient to work with dimensionless quantities. 
Rescaling the fields in (|I|) such that length is measured in units of the magnetic length 
Ib = I/a/c-B, chemical potential in units of the cyclotron energy ujc = 2eKB, time in units 
of the inverse cyclotron frequency and density in units of its constant bulk value, 

do = flo/^c, P' = fJ'/^c 

dt = dt/uj^ 

di = Isdi, di = Ibch 

P = p/p, 

one finds that the Lagrangian (rescaled by an overall factor 29/2ne^B'^) becomes 

1 



£ = 0* (z4 - ao - V^ + /i) - - I ( V - a)0 



-A 



' + ^e^^^S^^SA 



(7) 



with 



A 



At^e 



Thus, the theory contains only one free, dimensionless parameter (in addition to the chem- 
ical potential) |]14|. In particular it is worth noting that the statistics parameter 6 is 
absorbed in the rescaled interaction parameter A. 



3 Bulk properties 

We start by reviewing the main features of the CSGL model in the absence of an external 
potential, ^(x) = 0. It is known for the FQHE that at special fillings, v = l/{2m + 1), 
the ground state is characterized by a uniform electron density, p = ueB /2t:. As is well 
known [^, the CSGL model has this constant solution when the chemical potential takes 
a special value /io, 

(9) 
(10) 
(11) 

This solution simultaneously minimizes the kinetic term (-Dj0 = 0) and the potential 

n0)=/^i0p-Ai0r- 



_ Po 
P = P=2X = 


eB 

29 


a = — eA 




ao = 0. 





However, there exist other static solutions with the same constant density p away from 
the potential minimum. These solutions have oq 7^ and /i 7^ /io with /i — oq = /io- This 
is because the chemical potential enters the equations of motion only in the combination 
/i — ao; thus, a shift in the chemical potential can always be absorbed as a constant shift 
in oq. In this sense, the solutions of the equations of motion are independent of the value 
of /i. However, the energy does depend on yU as can be seen from Eq.(|^). This observation 
is useful for understanding how the QH plateaux emerge in this model, as we will discuss 
below (see also Curnoe and Weiss |1T5| ). 

The low energy charged excitations of the QH system are the localized quasiparticles 
and quasiholes. These are believed to be formed when the density is changed relative 
to the special filling p = p, e.g. by changing the strength of the magnetic field B. In 
the CSGL model these excitations are represented by vortex and anti-vortex solutions, 
whose existence has been shown in [|[ ^ and their properties examined by several authors 
T^, O. (We refer to the solutions with charge deficit as the vortex solutions and 



the solutions with excess charge as the anti- vortex solutions.) A general vortex solution 
is characterized by a topological quantum number s, which is s = ±1 for the elementary 
vortex/anti- vortex solutions. The corresponding charges are ±i/e, and the energies (relative 
to the ground state) we denote by e'" for the vortex and e"""" for the anti-vortex. 

For given s and /2, the dimensionless vortex/anti- vortex energies are functions of A only. 



These functions were studied numerically by Tafelmayer in Ref. [^, for s = ±1,±2 and 
p = Pq. As demonstrated by these functions there are two separate regions of A where the 
vortex solutions show qualitatively different behaviour. The regions are separated by the 
special point A = 1/2, referred to as the self-dual point. At this value of A the vortices 
are solutions of a linear (self-dual) differential equation, and the vortex energy is exactly 
proportional to the vorticity. Thus, the vortices are non- interacting. When A > 1/2, 
the energy of n vortices with topological charge 1 is lower than that of one vortex with 
vorticity n (as demonstrated in Ref. |I^ for n = 2). This means that a stable multi- vortex 



configuration with s = 1 can exist, since the vortices repel each other. For A < 1/2 the 
situation is opposite, so vortices attract each other; n vortices with s = 1 would tend to 
combine to one large vortex. We may then identify the interval A > 1/2 as physically 
relevant for the FQH effect, since only in this interval a stable density of vortices will be 
formed when the filling fraction is decreased relative to the plateau value l/{2m + 1). For 
anti- vortices there is no self-dual point. They have a repulsive interaction for all values of 
A. 

The existence of finite-energy vortex/anti- vortex solutions together with the freedom to 
vary the chemical potential for a given solution leads to the incompressibility of the ground 
state and the existence of QH plateaux within the CSGL model. We argued above that 
a given solution of the equations of motion is independent of the value of p, whereas the 
energy of the solution does depend on p. The energy (^ can be rewritten in a dimensionless 
form as 

E{p) =Eo- pQ, (12) 



where Q = 26Q = 26 J d'^x p, fl = fi/2eK,B, E = 26E/2eK,B, and Eq is independent 
of fi. (The occurrence of the additional (dimensionless) factor 26 in the rescahng of the 
total energy and charge is due to our choice to include this factor in the rescaling of the 
density.) For fl = fi^ the energy is minimized by the constant solution p = 1, and the 
vortex- and anti- vortex solutions have positive energies, e"{jlo) and e"''(/io) respectively, 
relative to the ground state. When the chemical potential is increased by 5/i, the energies 
of all the solutions will decrease by 5ji times their respective charge. Since the anti-vortex 
configuration has an excess charge as compared to p = 1, its energy will be lower than that 
of the uniform solution when fi exceeds a critical value 

/i+ = /io + e""(/io)/2vr (13) 

with 27r as the (rescaled) charge of the anti- vortex. Thus, for /i > /i+ the uniform ground 
state becomes unstable with respect to anti-vortex formation. Similarly, a lower critical 
value exists, 

/i- = /io-r(/io)/27r, (14) 

and if this is exceeded the uniform ground state becomes unstable with respect to vortex 
formation. Thus, there is a total "window" in the (dimensionless) chemical potential, 

A/i=(e-(/io)+e'^'^(/io))/27r (15) 

within which the ground state is given by the uniform density solution. 

Since fi = n/2eKB, a variation in fi can be interpreted either as a variation in /i at fixed 

Sp = ^5p (16) 



or as a variation in B at fixed /x. 



* = -2^^^- (^^' 



corresponding to the following interpretations (returning to dimensionful quantities): 
1) The ground state density at fixed B remains unchanged within a finite interval 

A/i=(e^(/io) + e"'^(/io))/^^ (18) 

which implies the incompressibility condition 

= (19) 




near /i = /iq. 

2) The density remains "locked" to the magnetic field through the condition b + eB = 

when B is varied within a finite interval 

AB = (e'^l/xo) + e-(/Xo)) -\ = (e-(^o) + ^"''(^o)) 7^. (20) 

ez/^A 47rA 

Since the Hall conductivity is proportional to p/B this implies that an is constant within 
the interval AB, giving rise to the Hall plateaux. 

So, at least qualitatively, the simple mean field theory reproduces the right picture as 
far as the plateaux are concerned. This is the case even without an impurity potential, 
since this is a model with a fixed chemical potential rather than fixed charge density. Thus, 
on a plateau the charge density changes with the magnetic field until it is energetically 
favoured to create vortices or anti- vortices in order to re-adjust the mean charge density. 
At this point the Hall conductivity leaves the plateau. 

When the strength of A is increased, the width of the plateaux will decrease. This 
follows from the expresssion (J2^) and the fact that the dimensionless vortex-and anti- 
vortex energies divided by A decrease with A (see Ref. [[l^). This is in accordance with 
the expectation that an increase in the (Coulomb) interaction of the particles will tend to 
suppress variations in the charge density. If an impurity potential is added we expect a 
broadening of the plateaux beyond the point where vortices are formed, since the vortices 
will bind to the impurities and therefore not affect the Hall conductance. This gives the 
connection to the standard picture of the fractional quantum Hall effect. 

4 The edge 

In the presence of a confining potential V^(x), corresponding to an external electric field 
eE = —W, the size of the system becomes finite, with a density profile at the edge that 
depends on the form of V. We choose here a geometry such that V^ is a function of x only 
and thus translationally invariant in the y-direction. This means that the edge is parallel 
to the y-axis. We further choose x = to be an interior line in the bulk (far away from 
the edge), and assume the system to be symmetric about this line. Due to this symmetry 
it is sufficient to consider the system for x > 0. We first examine the ground state of the 



system and then, in Sec. ^72|, discuss the gapless edge modes obtained by perturbing the 
ground state. Static solutions have previously been obtained numerically by Orgad and 
Levit IP in the case of an infinitely high step function potential. In this paper, we consider 
a smooth potential V{x) which is more relevant for real physical situations. Q 

^However, the translationally invariant form we assume for the ground state density at the edge may 
implicitly mean that the steepness of the potential should not be too small. A slowly varying potential may 
be considered as an x-dependent chemical potential, and from the discussion of the bulk properties, we in 
this case expect a broad edge with vortices trapped in the edge region. This corresponds in the FQHE to 
a composite edge with strips of compressible and incompressible Hall fluids, as discussed in Ref. p^. 



4.1 Ground state 

We consider ground states which are translationally invariant in the y-direction, i.e. they 
correspond to static solutions of the field equations where all the fields are functions of x 
only. The field equations then reduce to the form, 

Kdl(l){x) = (ao{x) - fi + V{x) + KAy{xf + 2Xp^(f){x) (21) 

d,Ay{x) = -2e{p{x)-p) (22) 

(9^ao(a;) = AK9Ay{x)p{x), (23) 

where we have made the gauge choice = y^, i.e. the wave functions are real. Again, 
Ay = cAy + ay and due to the symmetry of the system (and our gauge choice) a^, = 0. 
As before, A^. = and Ay = Bx. We further choose aj/(0) = 0, and as follows from 
the discussion in the previous section ao(0) = jj, — Hq va. order for the density to take the 
constant value p in the bulk. 

There is a conserved charge in this problem, Q = J dxp, and consequently a one- 
parameter set of ground states parametrized by this charge p . This set of ground states 
can alternatively be parametrized by the chemical potential. Thus, the ground state energy, 
which has the /i-dependence E = Eq{Q) — fiQ, may initially be considered as depending on 
two independent parameters Q and p,. However, for the true ground state, which is found 
by minimizing the energy with respect to Q, 

mj^ . 0, (24) 

the chemical potential will be a function of Q, 

dEn 

From the discussion of the previous section we know that there is no change in the bulk 
charge for small variations in p, around /iq. This means that when the chemical potential is 
changed within the limits determined by the vortex- and anti-vortex energies, all changes 
in the ground state take place at the edge. In particular, the additional charge which is 
introduced by a change in the chemical potential, is confined to the edge. 

The energy is stationary with respect to variations in the wave function about the 
ground state. This is true, due to (|2^) , also for variations where the total charge is changed. 
For an infinitesimal translation 5x of the edge, with 5Q = p6x and Sp{x) = —p'{x)6x, one 
finds 

SEo = Sx (Xp'^ - f dx V{x)p'{x)] (26) 



which leads to 



^ = i^ + }L I dx V'{x)p{x). (27) 

2 p 



The same expression, but with the additional term ao(oo) on the right hand side can be 
derived from the equation of motion ( pi]) by multiplying it with (f)'{x) and integrating over 
X. This leads to the conclusion ao(oo) = 0. An interesting consequence of this is that the 
total current in the system is zero when fi = /io: Integrating Eq.(^) over x gives 



ao(oo) -ao(0) 



-29J, 



(2J 



where J is the total current (in the ^/-direction). Since ao(0) = /i — yUg we have 

/i - /io = 26 J. (29) 

Thus, the integrated current, which in general will have two contributions of opposite sign 
at the edge, will have an overall sign which changes at /i = /io- 
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Figure 1: Ground state edge profiles for A = 1 and V{x) = \ 9{x — b) ■ {x — 5)^. The curves correspond 
to (from above) the upper critical value /i+, ji — ^q and the lower critical value /x~. 

This change in sign of the current has an analogy in the microscopic picture of the real 
quantum Hall system: There are two contributions to the current. One is the drift current, 
caused by the drift of the cyclotron orbits in the external electric field. The other is the 
polarization current, due to the gradient of the density at the edge and with origin in the 
cyclotron motion of the electrons. These contributions will have opposite sign, and which 



9 



one will dominate is determined by the value of the chemical potential or, equivalently, by 
the charge collected at the edge. 

A numerical study of charge and current densities in the ground state is presented 
in Figs. ID and |^. In Fig. ^ the edge profile 0(x) is shown for three different values of the 
chemical potential in the case of a harmonic confining potential ^(x) = 0.5 6'(a; — 5)-(x — 5)^ 
where 6{x) is the step function, and A = X/{Ak6) = 1. The upper curve corresponds to 
fi = 3.5 Uc which is the upper critical value as estimated from the numerical values in 
Ref. 



14 



Note the excess charge accumulated at the edge. For /x larger than this value the 
ground state is expected to be an anti-vortex configuration in the bulk rather than a state 
with uniform density. The middle curve corresponds to /i = /io = 2 cUc whereas /i = 1.3 cJc, 
the estimated lower critical value, for the lowest curve. For smaller values of the chemical 
potential the ground state is believed to reorganize into a density of vortices. 
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Figure 2: Current densities corresponding to the profiles in Fig.H. The total (static) current is negative 
at /i < /io, zero at /i = /io and positive for /i > /io- 

The current density profiles corresponding to these three solutions are shown in Fig. ^. 
Note the change in direction of the total static current: For fi = 1.3 Uc < /io the current is 
negative; for n = hq the two contributions cancel (as anticipated in Eq. (|29|) ), whereas the 
positive contribution dominates when // > /io as in the last plot. 
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4.2 Edge excitations 



The linearized field equations for fiuctuations of the fields about the ground state have the 
form 



V X 5ai 

djao + dfSa^ 

dySao + dtday 

dtSp 



+ 



[6ao + dtSe + 2KAy {6ay - dy69) + 2X6p] y/p 
ao + V - iJ, + kA^ + 2Xp b^ 
-2e5p 

AkO [Ay5p + p {5ay - dyS6)] 
-AnOp {5a^ - dJO) 
2kV ■ [A6p + p{6a- VdO)] . 



(30) 
(31) 
(32) 
(33) 
(34) 



For low-frequency fiuctuations propagating along the edge we write, as a simple ansatz, all 
the fields as an x-dependent part modulated by a slowly varying {y, t)-dependent harmonic 
part, 

J{qy~-ujt) 



6p{x,y,t) 
6e{x,y,t) 



6p{x) 
—irj{x) 



^i{qy-ujt) 



6a^{x,y,t) = 6a^{x) e'^'^^-^'\ (35) 

The field equations can then be reduced to purely x-dependent equations with the pre- 
scription dy -^ iq and dt — > —iu. We consider g as a small parameter and expand all the 
fields in powers of q, e.g. 

6p{x) = ,5p(°) (x) + q 5p(^) (x) + q^ 6p^^^ (x) + ■ ■ ■ . (36) 

Since we are interested in the gapless modes we assume uj to be of first order in q. The 
field equations (|5nD-(piD are examined order by order in q. As a starting point we assume 
that to lowest order only the phase fiuctuation rj contributes. This reduces the O*'* order 
equations to 

AKeipd^T]) = d, ipd^T]) = (37) 

so 1]^^' is a constant. 

only the terms involving p{5ax — d^O) oc Sj^ survive. 



<5j(i) = 0. (38) 

on the other hand, take the same form as those for 

3|) for the Q-dependent ground 



To first order in g, in (^3]) and 
simply giving 



The first three equations (|30|)-(|32| 

a first order variation in the equations of motion (|2T 

states. Thus, to 0{q) we have the following solution: 



6p^^'[x) 
6% (x) 

<5a«(x) 



[X 



6ai'\x)-6p+-v^^^ 



6al'\x)+r] 



(0) 



(39) 
(40) 

(41) 
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where ^^^ refers to variations between stationary fields obtained by a small change SQ in 
the ground state charge and a corresponding change 5/i in the chemical potential. 

For the stationary ground state fields the current conservation is trivially satisfied. 
This is not the case for the modulated fields (|35|), where the current conservation in fact 
determines the dispersion law for the edge waves. We note that to lowest order the current 
conservation involves the second order contribution to j^ and can be used to determine 
this from the first order fields. However, the integrated equation contains only first order 
fields, due to the boundary condition 

J.(0)=j.(oo) = (42) 

and has the form 

uj5Q = qSJ. (43) 

This equation is in fact correct to all orders in q, with 6Q as the fiuctuation in the integrated 
charge and 6 J as the fiuctuation of the total current in the ^/-direction. To first order in q 
this gives 

u dJ 1 d/i , , 

g ~ dQ~2edQ ^ ' 

where the derivatives refer to the Q-dependent ground state. This expression for the 
dispersion law agrees with the expression found in Ref. in the case of an infinitely steep 
confining potential. In our case we can further express the dispersion in terms of the 
potential V{x). Inserting the expression (^) for /i derived in the previous section into ( ^41) 
gives 

- = — dxV'ix)-^. (45) 

q en Jo dQ 

For a confining potential which is approximately linear over the edge it has the form 

with eE = —V, and more generally it can be read as an averaged form of this equation, 
with the ii^-field averaged over the edge. We note that the velocity of the edge waves is 
always in the same direction as the drift velocity of the electrons, and is not necessarily in 
the same direction as the (integrated) edge current in the ground state. 

5 Two-component system 
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5.1 The model 



The CSGL Lagrangian (P can be extended to include spin |[T^. This is done by replacing 
by a two-component complex scalar field and introducing a Zeeman term 



'~'Zeem „ *?* '^ z 



(47) 



The model then takes the form 

2 



-^ = H <^a ^^^i -a^^ y^(x-y) 



a=\ 






A(pi + 92? + —e^'^^a^d.ax, 



(48) 



where the index a denotes the two spin states (up and down), so that p = pi + P2 is the 
total (charge) density, while ps = Pi — P2 is the spin density. A spin unit vector may be 
introduced by 



n 



1 

P 



'crc 



(49) 



with components in the plane as well as in the direction of the -B-field (z-direction). In 
(^) the Zeeman term is represented by the difference between the chemical potentials of 
the two components. 

With the densities and phases of the two components of (pa as variables, the field 
equations are. 



«V2^/p^ =\aa + V- (p, - dtda) + k{A- V9af + 2A (pi + pa) 

V X a = -2^ (pi + P2) 

•9x00 + dta^ = Ak9 J2a Pa {Ay - dy9a) 

dyao + dttty = -Ak6 J2a Pa {ax - d^Oa) 

dtpa + V ■ j„ = 0, 

where Ai = ai + eAi as before, and the current densities are given by 

j„ = -2Kpa {A - Vdo) . 

This system has two conserved charges Qa = / (Px pa ■ 



'Pa 



(50) 



(51) 



5.2 Ground state properties 

The form of the Hamiltonian 

2 

'^ = >^^\{^ - iA)(l)a\^ + V{pi + P2) + A(pi + P2f - PlPl - P2P2 

0=1 



(52) 
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implies that the bulk ground state {V = 0) is fully polarized and has a uniform density. 

^eB _ Ho 

0, 



Pi 

P2 



^ 29 2X 



(53) 
(54) 



where we have assumed fii > fi2- The chemical potential of the upper component is related 
to iiQ by 



Hi - ao(0) = Ho- 



(55) 



Even though the ground state in the bulk is fully polarized, this may not be the case at the 
edge of the system. Tilting the spins {i.e. p2 ^ 0) at the edge may lower the energy if the 
Zeeman gap is not too large. This problem was recently studied numerically in Ref. |I3[ ; 
their calculations were based on Hartree-Fock and an effective action approach |T^ in the 
low-energy limit of the CSGL model. Here, we shall examine some of the ground state 
properties at the edge both analytically and numerically, within the framework of CSGL 
theory. 

As before, we introduce a confining potential V{x), which gives rise to a straight edge 
parallel to the y-axis. We consider stationary solutions, with densities pi and p2 that are 
translationally invariant in the y-direction. With the gauge choice 6i = 0, the solutions 
take the general form 



pi(x) 



p2[x)e 



i{ky-f3t) 



(56) 
(57) 



with pi = p, p2 = in the bulk. The gauge choice implies A^ = a^^ = with oq and ay 
as functions of x only. The parameters k and (3 are undetermined at this stage. However, 
we see that if fc 7^ this means that the spins will rotate around the i?-field when moving 
along the edge. With the ansatz ( |56D and (|^) the equations of motion reduce to 



i^dl^2 
dxtto 



ao + V - pi + kAI + 2A(pi + P2)] ^/p^ 

ao + V - p2- P + n{Ay - kf + 2A(pi + pa) 
-26 (pi + P2 - p) 
AKe[piAy + p2{Ay-k)] 



(58) 

(59) 
(60) 
(61) 



The corresponding energy (per unit length in the y-direction) is 



E 



dx 



1^ {dxy/pif + /t {dx^/p2f + k^AIpi + K{Ay - kfp2 + V{pi + P2) 



A(pi + P2)^ -piQi -P2Q2, 



(62) 
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where Qa now denotes charge per unit length, Qa = J dx Paix). 
The ground state conditions are 



dE dE 

(63) 



dQi dQ2 

and 

These conditions determine the chemical potentials as functions of the charges, fia = 
fJ'a{Qi,Q2)- So here we have a two-parameter set of ground state solutions parametrized 
by the charges Qi and Q2. 

From the condition (|63D, that the ground state energy is invariant under any first order 
variation of the two charges, we can derive expressions for pi and f3. We first consider 
a perturbation in Qi caused by an infinitesimal translation of both components in the 
x-direction (thus keeping Q2 fixed). This gives 

ui = Ap + - dx V'(x)p(x) (65) 

pJ 

which corresponds to Eq.(^) in the one-component case. Again, the same expression but 
with an additional term ao(c)o) can be derived directly from the equations of motion, giving 
ao(oo) = in the ground state, and consequently 

Pi — Po = 29 J, (66) 

where J = Ji + J2 is the total (integrated) current along the edge. 

Similarly, we can slightly perturb the "spin charge" Qs = Qi — Q2 keeping the total 
charge fixed. This is done by an infinitesimal SU{2) rotation of (p, 

01 — ^ 01 + e4>2 (67) 

02 -^ 02-£*0i (68) 

with 

e = e{y) = e,e-'^y (69) 

and £0 as an infinitesimally small real constant. This introduces the following variations 
in the charges. 



6Qs = 4:eoJ dx^pi{x)^p2{x), SQ = 0, (70) 

and the condition of stationary energy, dE/dQs = 0, gives 



Jdx ^pi{x)^p2{x) 
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The same expression, but with an additional factor (3 on the right hand side, can be deduced 
from the equations of motion by multiplying (^) and (^9]) by ^/p2 and ^/pi respectively, 
subtracting the two equations and integrating over x. This leads us to the conclusion that 
the ground state condition dE/dQs = implies /3 = 0. Thus, there is no t-dependence 
in the phase of 02- Note that the relations (|66|) and (|7l|) were derived only demanding 
that the energy be minimized with respect to the charges; they are therefore valid even for 
values of k that do not minimize the energy. Finally, also demanding dE/dk = one finds 



Ja = -2k / dx p2{x) {Ay{x) -k) = (72) 



which gives an expression for the ground state value of k, 



k = ko = -—dx p2ix)Ay{x). (73) 

As already stated, we expect a critical value of the Zeeman gap A/i = pi — fi2 to exist, 
so that if this critical value is exceeded the system will be fully polarized also at the edge. 
If the gap is close to, but smaller than this value we expect the second component of the 
wave function to be small. We will make an analysis of this case and illustrate the case by 
numerical evaluations of the ground state wave function. 

For small Q2 we can expand p2 in powers of this charge (for given Qi and k) 

/X2 = /if +4'^Q2 + --- (74) 

and similarly for the density p2, 

P2 = P2/Q2 = P^?+P?Q2---. (75) 

The lowest order terms correspond to the limit (^2 — ^ 0. This means that p2 is the value of 
P2 which corresponds to the critical Zeeman gap for given pi. Similarly, pg represents the 
profile oi P2 in the same limit. In this limit p2 is negligible compared to pi, and the problem 
is simplified since the equations for the two densities ((|58|)-(pTD) then decouple: Eqns.(|5^), 
(|60|) and (|6TD get reduced to the single-component equations (pTD-(p3D discussed in Sec. ^J] , 
and the equation (EDI) for ^/p^ simply becomes a one-dimensional potential problem. 



'^9l^p2ix) = [U{x, k) - P2] ^p2ix) (76) 

where the potential is 

U{x, k) = ao{x) + V{x) + K{Ayix) - kf + 2Api(a;). (77) 

The ground state is found by minimizing the energy with respect to k, and the chem- 
ical potential p2 is determined as the (lowest) eigenvalue of the potential problem. This 
determines the critical value of the Zeeman gap for the chosen value of pi. 
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Figure 3: The (dimensionless) potential U{x, fco) determining the density p2{x) of the lower component 
in the limit p2 ^ pi- Here, the confining potential is V{x) = 0.05 9{x — 5) ■ {x ~ 5)^, pi w 1.5 and 
ko = ko-lB = 0.49. 



The expansion in Q2 may be used to correct iteratively the lowest order results. Thus, 
the lowest order contribution to p2 gives first order corrections to the equations for pi, Oq, Ay. 
This in turn gives first order corrections to the potential U{x, k) and to fi2, P2 and ko, etc. 
This gives a systematic way to determine the ground state density profiles for small Q2- 

We have solved this problem numerically only to lowest order in Q2. Harmonic poten- 
tials of varying strength have been used, and several values for /xi have been chosen in the 
allowed interval (as determined in the one-component case). 

Fig.|^ shows the potential U{x,ko) for one specific value of /ii. The minimum of U is 
located at the edge as expected. The wave number k^ takes a non-vanishing value. Thus, 
the spin vector is rotating when moving along the edge. 

The corresponding density profiles are shown in Fig. ^ The figure shows that the lower 
component {ip2) of ground state density is relatively weakly bound in the potential U{x, k) 
for this particular choice of V{x) and /ii. The localization gets sharper for values of /ii 
corresponding to larger k^. A similar effect is seen in Fig.|^ when k is changed relative to 
the ground state value ko. 

Table |] shows numerically estimated values of ko and the critical Zeeman gap for dif- 
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Figure 4: Edge profiles ^/pi(x) = (piix) and yjp2{x) = ^2{x) corresponding to the potential in Fig.|^. 
The corresponding (critical) Zeeman gap is /ii — /i2 ~ 0.002. 

ferent jii close to the lower critical value /i~ (discussed in Sec.|]), and for a fixed confining 
potential V{x). When /xi is increased towards /io, both the wave number and the critical 
Zeeman gap decrease until A/z^r goes to zero with a finite value of /cq. A corresponding 
study with a steeper confining potential shows no solution with positive k^ and A/ic,- close 
to /i~. This indicates that the interval of yUi where spin textures occur, increases with 
softening of the confining potential. 

Similarly, one finds solutions with A/Xcr > near the upper critical value /i"*". These so- 
lutions are characterized by negative wave numbers whose absolute value decreases along 
with the critical Zeeman energy as /ii is lowered. 

We can summarize our results as follows: For a sufficiently soft confining potential and 
for Zeeman energies smaller than some critical value, the spins at the edge are tilted in the 
ground state (n^ 7^ 1) and rotate around the z-axis with a wave number fco as one moves 
along the edge. When ni is increased from its lowest allowed value /i~ at fixed confining 
potential, both k^ and the critical Zeeman gap A/Xcr decrease until A/ic^ goes to zero. Near 
the highest allowed value /i"*" the edge is again spin textured, with negative /cq and a critical 
Zeeman gap which approaches zero as /xi is decreased. The ranges in /ii where the edge is 
spin textured, decrease as the confining potential gets steeper. 
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Figure 5: Shape of 1^92(2;) = \/p2{x) for values of k away from fco- The confining potential is V{x) — 
0.05 9{x — 5) • (a; — 5)^ and jli — 1.5. The curves correspond to (from below) k = 0.49 (~ fco); 0.55, 0.60. 



As mentioned in Sec. ^, the egde is expected to become unstable with respect to a redis- 
tribution of the charge at soft confining potentials; for /ii away from fiQ this reconstruction 
is expected to take place at steeper V{x) than is the case for /ii = /xq- This means that 
for fixed fii texturing may occur only in some intermediate region where the confining 
potential is neither too steep nor too soft. 

Since the "window" (Eq. (|T8|) ) of allowed /^i depends on the interaction strength A, 
the results discussed above may change quantitatively for other values of A; here we have 
only examined the case A = 1. For this value of A we have found values for the critical 
Zeeman gap (see table |l|) which are of the same order of magnitude as the Zeeman gap 
reported in real experimental situations; typically, Afi/uJc ~ 1/60 in QH experiments in 
GaAs heterojunctions [BO]. 



5.3 Edge excitations 

When the edge is spin textured one may expect new massless edge modes to appear; as 
argued in Ref. [0|, translational symmetry along the edge is spontaneously broken by the 
phase of 02, and one may expect a corresponding Goldstone mode to exist (at least at the 
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h 


ko 


After 


1.3 = jl- 


1.0 


0.06 


1.4 


0.6 


0.01 


1.5 


0.5 


0.002 


1.6 


0.3 


-0.007 



Table 1: Wave number and critical Zeeman gap A/i^. = /ii — /i2 for some values of the 
chemical potential /ii at A = 1. The confining potential is V{x) = 0.05 ■ 6{x — 5) ■ 
{x — 5)^. Chemical potentials are measured in units of the cyclotron energy, whereas the 
dimensionless wave number is k = k ■ Ib- 



mean field level). We have examined the edge modes by considering quadratic fiuctuations 
in the fields around the ground state configuration in the same way as was done for the 
spinless case. 

Again, we assume the fiuctuations in the fields to be separable in the x- and y- 
dependence, 



6p^ix,y,t) = 5p,(x)e^('?^--*) 
6e^{x,y,t) = -tr]^{x)e'^'^y-^'^ 
6a^ix,y,t) = 5a^(x)e^('?^--*) 

and thus replacing dy —>■ iq and dt —>■ —iuj, one finds for the linearized field equations 
K [dl - q^j 5^fp^ = [6ao - ujT]a + 2k {Ay - k^) {6ay - qr]a) + 2X6p] y/p^ 



(78) 
(79) 
(80) 



dx6ay — iqdttx 
dxSao — iuuSttx 

iqSao — iuoday 

—iujdpa 



+ aQ + V - pa + i^{Ay - kaf + 2Ap 6^0 
= -2e6p 

= Ak9 Y^ [5pa {Ay - ka) + Pa {^O-y " qVa)] 
a 

= -AkO J2 [paSa^ + ipadxVa] 
a 

= 2k [iq6pa {Ay - ka) + iqpa {^o-y - qVa) 
+ dx {pa {Sax + idxr]a))] 



(81) 
(82) 
(83) 

(84) 
(85) 



with ki = and k2 = k. As in Eq. (|36|) we expand all the x-dependent functions in powers 
of q and assume that the phases r/i and 772 dominate to O*'^ order. In analogy with Eq. (p7|) , 
one finds 



T.Pc^9xvl^^ = 9.{padxvl^'>)=0 



(86) 



which implies that r][ and 772 ^ire constants. 
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To 0{q), the last two equations ( |51D and ( P3| ) give 

s£}ot = dMl = 0, (87) 

and therefore Sj^^l^ = 0, as follows from the boundary conditions on the current deep in 
the bulk and outside the edge. 

The first three equations, (|8lD -(p3D, again take the form of those for variations in the t- 
and ?/-independent (ground state) fields {5p^\ ^o,}f\ ^l^a, Sk) (see Eqs. (p^)-(pT|)), provided 
we make the following identifications 



H^^-^r = 


- H'' 




H^^-^f = 


- H^^- 


-5k 




= Sai'^ 


-6 1^1 


X (1) ^ (0) 

Say --V2 -- 
q 


- Sai'' 


- 5^2- 



(88) 



One should note that in these expressions 6k is not restricted by the condition (0), only 
the ground state conditions (|63D are satisfied. Thus the variations in the fields depend on 
three parameters, which we may choose to be Qi,Q2 and k. However, the identifications 
(^) implicitly give an additional constraint, 

u!6k = —q{6fii — 6fi2) = —qSAfi (89) 

and with this condition included, the fields depend only on two parameters, Qi and Q2- 
The continuity equations for the currents, 

uj5Qa = qSJa (90) 

give two linear, homogenous equations in the two variables which determine the (low 
frequency) normal modes. The equations have the form 

M^^6Qf, = (91) 

with 

Mal3 = (- + Aflk)i-Saf3 - ^q/s) + Jak^fi^, (92) 

where we have introduced the notations 

^^' = ^k- ^^" = 9q: 

J.k = ^ U = ^. (93) 
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Eq.(PTD corresponds to Eq.(^^ for the one-component case, which there gave uj/q in terms 
of derivatives of J and /i with respect to Q. Here the matrix equation gives a secular 
equation which is third order in uj/q, 



' \ 3 / \ 2 

Uj\ . Uj\ „UJ 



- +A\-\ +5- + C = (94) 



with 



A = A/ifc - Jii - J22 

B = Jii J22 - ^12^21 - Ajj,k{Jn + J22) + A/ii Jife + Ajj,2J2k 

C = A.fik{JiiJ22 — J12J21) — A^i{JikJ22 — J2kJl2) — 'AfJ,2iJ2kJll — JlkJ2l)- (95) 

The parameters of this equation are determined by the solutions of the (ground state) 
equations (|6T]). We do not have any explicit expressions for the parameters, but one 
should note that they are not all independent, due to relations like (|66D . 

The fact that ( P^ is third order in u/q indicates that with spin effects included there 
are three massless chiral edge modes (rather than two). We cannot give the velocities of 
these modes in explicit form, but we note that the expressions (^) and (§§), which relate 
the total current to the confining potential, imply that the velocities can be written as 

where 6p and 6Q are the variations in total charge density and integrated charge, respec- 
tively, for each of the normal modes (provided 6Q ^ 0). When the variation in the charge 
density is located at the edge this shows that the velocities of all the modes can be written 
in the general form — < E > /B, with < £" > as an averaged value oi E = —V'/e at the 
edge. However, the ii^-field will be averaged with functions which in general are different 
for each of the three modes. 



6 Discussion 

The Chern-Simons Ginzburg-Landau theory gives a mean field description of the fractional 
quantum Hall effect which reproduces correctly several of its main features. Here we have 
examined, in particular, edge effects in this framework, both by analytical and numerical 
methods. The starting point has been the observation that if the chemical potential is 
varied within a certain interval there is no change in the uniform bulk density, but a 
change in the edge profile is introduced. Thus, there is a set of ground states parametrized 
by the edge charge. The edge current also varies with the chemical potential, and we have 
found that it changes sign for the value of the chemical which minimizes the potential 
energy in bulk. This value can be considered as corresponding to the central point of the 
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Hall plateau, the point where the uniform charge density of the electrons neutralizes the 
(constant) background charge. 

The low-frequency edge modes have been examined in terms of quadratic fluctuations 
in the fields and with an expansion in the wave number q for fluctuations propagating along 
the edge. As expected (and previously shown for an infinitely steep confining potential) 
there exists in the spinless case one gapless chiral edge mode. To first order in q the density 
profile of this mode is the same as for an (infinitesimal) variation in the ground state profile 
introduced by a variation in the charge (and the chemical potential). An expression for 
the velocity of propagation has been found in terms of the confining potential; it has the 
form of an averaged drift velocity E/ B^ averaged over the edge with the density profile of 
the gapless mode. 

The edge effects of a fully polarized system with spin have been examined in a similar 
way. In this case there will be a spin texturing of the edge if the confining potential is 
sufficiently soft and if the Zeeman energy is not too large. The ground state then depends 
on two edge charges, one for each spin component, and we have examined the case were the 
second charge (corresponding to the spin down component) is small. The critical Zeeman 
gap has been examined numerically as a function of the chemical potential (for a particular 
value of the interaction strength). It is different from zero in intervals of ^i near the upper 
and lower bounds /i"*" and /i~ and decreases when /xi approaches /iq from either side. The 
width of these intervals increases with softening of the confining potential. 

Also in this case there are gapless edge modes which can be related to variations in the 
ground state fields introduced by variations in the charges. However, the variations involve 
changes in the wave number k of the spin down component which are not restricted by the 
ground state condition dE/dk = 0. 

The continuity equations for the charges give a third order equation for the mode 
frequencies. This indicates the presence of three chiral modes. However, the equation has 
been given in terms of a set of undetermined parameters. It will be of interest to examine 
further these parameters to see if there are more relations or constraints than have been 
established in this paper. Also a numerical study of these modes is an interesting subject 
for further research. 

The edge modes have here been studied as quadratic fluctuations about the mean field 
ground state. A question which deserves further study is whether quantum effects will 
qualitatively change this picture. The reason that this cannot be ruled out is the quasi- 
one dimensionality of the system and the suggested connection to spontanous breaking of 
translational symmetry of the ground state [^. As is well known, spontaneous symmetry 



breaking is not supported in one-dimensional systems due to (quantum) fluctuations in the 



symmetry breaking variable [22 
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